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Abstract 

We start by presenting physical arguments for the impossibility of perturbative ther- 
malization leading to (non-viscous) Bjorken hydrodynamic description of heavy ion col- 
lisions. These arguments are complimentary to our more formal argument presented in 
PP. We argue that the success of hydrodynamic models in describing the quark-gluon 
system produced in heavy ion collisions could only be due to non-perturbative strong 
coupling effects. We continue by studying non-perturbative effects in heavy ion collisions 
at high energies. We model non-perturbative phenomena by an instanton ensemble. We 
show that non-perturbative instanton vacuum fields may significantly contribute to jet 
quenching in nuclear collisions. At the same time, the instanton ensemble contribution 
to thermalization is likely to be rather weak, leading to non-perturbative thermalization 
time comparable to the time of hadronization. This example illustrates that jet quenching 
is not necessarily a signal of a thermalized medium. Indeed, since the instanton models 
do not capture all the effects of QCD vacuum (e.g. they do not account for confinement), 
there may be other non-perturbative effects facilitating thermalization of the system. 



*e-mail: yuri@mps.ohio-state.edu 



1 Introduction 



The problem of thermalization has become one of the central questions in the theory of heavy 
ion collisions. The success of hydrodynamic models [3 13 13 13 E] in describing RHIC data on 
elliptic flow [7j, along with the suppression of high-j9 T particles in Au+Au collisions as compared 
to d + Au and pp collisions |Hl 13 U01 HU [T3 H31 [HI |T3l HZl UHl HH1 UHl I20| , indicate creation 
of a strongly interacting and likely thermal medium in the final state of heavy ion collisions 
- the quark-gluon plasma (QGP). Currently, the main theoretical puzzle is in understanding 
how this medium was formed in the collision. 

The initial conditions for creation of QGP has received a lot of attention in the recent 
years [23 122 123 121 123 1213 123 in the framework of saturation/Color Glass Condensate (CGC) 
physics I23I23E3E3E3E3E11E3E3EI1EE1- The data on hadron multiplicity and spectra 
obtained in d + Au collisions at RHIC [23 113 EJ H3 ESj appears to be in good agreement with 
expectations of Color Glass approach HH EH1 HHl HZl EHl EHl EIIl EH E21 EB] When combined 
with the success of saturation models [HI] in describing total charged hadron multiplicity in 
Au + Au collisions at RHIC, this allows one to hope that evolution of the produced quark- 
gluon system toward thermalization could also be described within the saturation/Color Glass 
framework, which indeed involves using Feynman diagrams. 

However, understanding perturbative thermalization proved to be a very hard problem. 
Several perturbative thermalization scenarios have been proposed in the recent years [S3 E3 
I5T1 E3 E3 E3 EI] ■ None of them appears to lead to thermalization time as short as required 
by RHIC phenomenology, which is r<0.5 fm/c. Moreover, the so-called "bottom-up" thermal- 
ization scenario, which was originally proposed in [55] and was the first one to use saturation 
initial conditions, proved to be susceptible to plasma instabilities [FS ] IM | IfiK ] Ififi ] IH7 | IfiH] . 
This brought hope that instability effects would speed up the onset of thermalization. How- 
ever, in jnn] it was shown that instabilities (along with any other dynamical effects) obey a 
lower bound on thermalization time, which is parametrically very close to (though slightly 
lower than) the thermalization time obtained in the "bottom- up" scenario [S3]. More recent 
developments include suggestion to consider the onset of isotropization instead of thermaliza- 
tion as the necessary condition for hydrodynamic models to work [HE], a numerical analysis of 
instabilities demonstrating that their magnitudes saturate as time progresses [10], leading to, 
probably, a less severe impact on the system as thought previously, and a new modified version 
of the "bottom-up" scenario [7T] . 

In our recent paper on the subject we considered the onset of thermalization as a possi- 
ble transition from the e ~ 1/r scaling of the energy density characteristic of the Color Glass 
condensate initial conditions [23 [23 123 123 123 123 EI], to a steeper falloff with proper time, 
e ~ 1/t 1+a , expected if the system obeys ideal Bjorken hydrodynamics |2 j (A > 0). We were 
looking for perturbative corrections, which could be, for instance, powers of A lnr, with A 
being proportional to some power of the strong coupling constant a s . Having found no such (or 
any other) important corrections to 1/r scaling of energy density we proved that no corrections 
of this type are possible in perturbation theory. Therefore, we proved that perturbative thermal- 
ization is impossible in heavy ion collisions. Our arguments in ^ were somewhat involved, and 
lacked a clear physical argument demonstrating impossibility of perturbative thermalization. 
Below we try to construct such an argument in Section 2. There we provide two arguments. 
The first one, presented in Sect. 2.1, demonstrates that, in the absence of non-perturbative 
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effects, if the theory is defined only perturbatively, Bjorken hydrodynamics would not lead to 
decoupling of the system. Therefore, plane-wave final states necessary in perturbation theory 
to define scattering amplitudes and particle production cross sections are never achieved. This 
indicates a contradiction between Bjorken hydrodynamics and perturbation theory, indicating 
possible impossibility of obtaining the former from the latter. The discrepancy becomes appar- 
ent in the second argument, presented in Sect. 2.2. There we argue that, in the zero coupling 
limit (g — > 0) taken at some late time after the collision, QCD perturbation theory would 
indeed lead to "free-streaming" with the energy density scaling like e ~ 1/r, while Bjorken hy- 
drodynamics would give e ~ 1/r 4 / 3 . The formal reason for such scaling in the case of Bjorken 
hydrodynamics is that the (non-interacting) ideal gas equation of state is e = 3p, and this is all 
that's needed to obtain the e ~ 1/r 4 / 3 scaling. The physical reason for e ~ 1/r 4 / 3 scaling in the 
g — > limit is that, in writing e = 3p equation of state for the ideal gas, one indeed assumes 
a presence of some thermal bath (e.g. the box containing the particles) with which the gas 
particles interact and which keeps them thermal throughout the evolution of the system. Such 
external thermal bath is absent in heavy ion collisions. Therefore, Bjorken hydrodynamics is 
impossible to achieve in the weak coupling limit. We further argue that another way to reach 
this conclusion is by noticing that shear viscosity, which is the first non-equilibrium correction 
to Bjorken hydrodynamics, diverges in the weak coupling limit (see Eq. (|27|)) 1721. indicating 
that non-equilibrium corrections become comparable to the equilibrium values of the pressure, 
breaking down the ideal (non- viscous) hydrodynamics in the weak coupling limit. Our con- 
clusion here is that thermalization of the quark-gluon system (if it does take place) leading 
to hydrodynamic description of the heavy ion collisions could only be due to non-perturbative 
strong coupling effects. 

Indeed, understanding the non-perturbative phenomena in the non-equilibrium setting of 
heavy ion collisions is a very difficult (if not presently impossible) task. Here we will consider 
one well-studied model of non-perturbative effects — the instanton vacuum model [HU EU [7J3 
EH1 EZl EH EH] • Since, as we argue in Sect. 2, non-perturbative effects appear to be important in 
AA collisions for thermalization, which is required by phenomenology to describe the observed 
data on elliptic flow |H E] , we would like to investigate their influence on other important 
observables in heavy ion collisions. In Sect. 3 we study the effect of the instanton vacuum on jet 
quenching and thermalization in heavy ion collisions. Reviewing the results of [HHj in Sect. 3.1 
we argue that the instanton density is suppressed by strong perturbative gluon fields produced 
(in the Color Glass Condensate framework) in the early stages of the collision with the proper 
times t<1/Q s . At slightly later proper times, for t>1/Q s , perturbative fields linearize and 
instanton vacuum is not suppressed anymore. At this point we would have a system of free 
quarks and gluons propagating through the instanton vacuum. Usually, a parton produced 
in a collision and propagating as a free plane wave would start hadronizing right away. This 
is probably true in pp or e + e~ collisions. However, in AA collisions, at times t>1/Q s , the 
energy density is still too high for hadronization to take place. Unable to hadronize due to high 
energy density, non-interacting quarks and gluons are likely to interact with the instantons in 
the vacuum and to loose energy as they propagate through the instanton ensemble until the 
hadronization time. (Once quarks and gluons are converted into hadrons they, of course, do 
not loose energy in the vacuum anymore.) The mechanism of parton-instanton interactions has 
been extensively studied in |81j . One might expect the energy loss of partons in the instanton 
vacuum to be negligibly small due to low density of instantons, which would play the role of 
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scattering centers. However, the peculiar property of the energy loss formula derived in |H2*1 IHB"j 
and shown here in Eq. ()47|) is that the integrated partonic energy loss is particularly sensitive to 
late proper times of the collision due to an extra factor of r in the integrand. Namely, it is not 
as important what the energy density was in the early stages of the collision, as it matters what 
the energy density is at late proper times. Therefore, the energy loss of a fast moving parton 
in a non-static medium is mostly determined by the medium density at late times. Since, at 
late proper times close to the time of hadronization, the density of scattering sources in the 
produced medium is comparable to the density of instantons in the vacuum, the latter may 
significantly contribute to jet quenching, on par with traditional energy loss in the produced 
medium, as we demonstrate by parametric estimates in Sect. 3.2. 

At the same time, as the vacuum instanton density is low, its contribution to thermalization 
appears to be rather weak, leading to non-perturbative thermalization time comparable to the 
time of hadronization, as we show in Sect. 3.3. This example illustrates that jet quenching is 
not necessarily a signal of a thermalized medium and, if Eq. (|47jl adequately describes energy 
loss in a heavy ion collision, one may have significant jet quenching without thermalization. Of 
course, as we argue in Sect. 4, since the instanton models do not capture all the effects of QCD 
vacuum (e.g. they can not explain confinement), there may be other non-perturbative effects, 
such as QCD strings [HI], facilitating thermalization of the system and leading to energy loss 
in the non-perturbatively formed QGP. 

2 Non-perturbative Nature of Bjorken Hydrodynamics 

In this Section we will present two qualitative arguments showing that perturbative thermal- 
ization leading to Bjorken hydrodynamics description of heavy ion collisions is impossible. 

2.1 Argument I 

Before we begin, let us briefly review the main assumptions of Bjorken hydrodynamics 2\. The 
main idea behind the ideal (non-viscous) Bjorken hydrodynamics is that the medium produced 
in a collision of two very large nuclei can be described by hydro dynamic equations 
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with the energy-momentum tensor given by 
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where, in the boost-invariant (rapidity-independent) case the fluid velocity is 




u = 
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and the proper time is r 



^J2x^xZ. We use the light cone coordinates 



t±z 
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and u denotes a two-dimensional transverse vector, e is the energy density and p is the pressure 
of the fluid. Substituting Eq. (J2J into Eq. ((TJ) and assuming that e = e(r) and p = p{r) yields 

m 

de e + p 



dr r 

Eq. (JHJ) leads to the scaling of energy density with the proper time [2j 



(5) 



(6) 
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where 



r l+A 



< A < |. (7) 

The value of A = 1/3 can be obtained in the ideal (ultrarelativistic) gas case with the equation 
of state 

e = 3p (8) 



and leads to 



The value of A = 0, leading to 



e (10) 

T 

could in principle be achieved in a QCD plasma near the critical phase transition temperature 
T c if the phase transition was first order j2j |EE]. However, the confining phase transition in 
QCD is non-perturbative: therefore it is impossible to obtain QCD matter at T = T c through 
Feynman diagrams and we will not be concerned with this case here. 

Now, let us assume that, as a result of summation of perturbative diagrams to any order 
(and to all orders if necessary) one shows that in a heavy ion collision a thermalized medium is 
produced which is describable by Bjorken hydrodynamics. With logarithmic (due to running 
of the coupling) accuracy the energy density of the resulting medium would scale with the 
temperature T as 

e ~ T 4 . (11) 

This formula is better justified in the weak coupling limit, where the diagrammatic analysis is 
at its best too. Comparing Eq. (fTTf to Eq. (JHJ) yields 

T ~ -4a- (12) 

T 4 

For example, for ideal gas Eq. ([T2|) yields T ~ t -1 / 3 . 

Let us now imagine that thermalization and the subsequent evolution of the system are 
completely described by the perturbation theory. In other words, let us imagine "turning off" 
confinement and all other effects associated with the non-perturbative QCD scale Aqcd an d 
assume that the theory is defined perturbatively. Indeed this is not a realistic assumption, but 
would help us clarify whether thermalization in heavy ion collisions is of perturbative nature 
or not: if thermalization is perturbative we should not get any inconsistencies by "turning off" 
confinement effects. If we do get inconsistencies, then perturbative thermalization would be 
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impossible, and if thermalization does take place in heavy ion collisions, it could only be due 
to non-perturbative effects. 

For simplicity let us consider only gluons and massless quarks in the system. Since the 
temperature T is now the only scale in the problem, the perturbative mean free path of a 
(massless) quark or a gluon in the medium would scale like 

Xpert ~ (13) 

which, using Eq. (|T2*j) gives 

Xpert ~ r 1 ^. (14) 
Employing Eq. (|ZJ) we obtain an upper limit on the mean free path 

X P en ~ r 1 ^ < r 1 / 3 . (15) 

The equality in Eq. (|15|) is achieved only for an ideal gas with the equation of state given 
by ©• (In Eq- (jUl) we have omitted r-independent factors which would make the inequality 
dimensionally correct: putting them back into the formula would not modify the conclusion, 
but would unnecessarily complicate the equation. The same applies to Eq. (fTTj) below.) 

Noting that the typical longitudinal extent of the system at time t after a heavy ion collision 

is 

L = t > r (16) 

we conclude that 

Xpert < T 1/3 < T < L. (17) 

(The transverse extent of the system is assumed to be infinitely large [2] and is, therefore, 
always larger than the mean free path.) 

Eq. (|17|) implies that the mean free path of a gluon or a massless quark in a plasma de- 
scribed by Bjorken hydrodynamics is always smaller than the size of the system. Therefore, 
the decoupling condition 

A ~ L (18) 

is never achieved in perturbative Bjorken plasma. 

(Indeed, once we turn the confinement effects back on, we would get another scale in the 
problem, Aq C d, which would lead to 

\ 1 1 3(1+*) / 1n x 

Xnon-pert = ~ 1 ~ T 4 (19) 

1V QCD 

giving the mean free path 

Xnon-pert ~ T ~ L (20) 

for the ideal gas case (A = 1/3), which marginally satisfies the decoupling condition. Here a 
is the parton-parton scattering cross section and p is the density of matter.) 

As we can see from Eq. (|17p. the asymptotic r — > +00 state of a perturbative plasma de- 
scribed by Bjorken hydrodynamics is not "free streaming": in fact, the inequality in Eq. (JT7j) 
gets more justified as r increases. Each individual parton can escape the system with an ex- 
ponentially small probability, e _L / A ~ exp(— const r 2 / 3 ), which only decreases with time. This 
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behavior is in contradiction with the diagrammatic approach: the final state of any Feynman 
diagram in high energy scattering is a set of free non-interacting particles flying off to infinity. 
The existence of such a plane wave final state is an essential assumption of the perturbation 
theory. Without the plane wave final state the problem of particle production would be ill- 
defined. If diagrams would lead to Bjorken hydrodynamics, than, such a common observable 
as the one-gluon inclusive production cross section would be impossible to define and calcu- 
late in AA collisions, since there will be no plane-wave particles in the final state. Single 
gluon inclusive production cross section is well-defined in the quasi-classical approximation of 
McLerran-Venugopalan model [HH EH 1221 EH EH EE] , which does not led to a thermalized final 
state |24j . If thermalization comes in through the higher order quantum corrections to classi- 
cal fields jnSj, such corrections would invalidate the notion of one-gluon inclusive production 
cross section, which appears absurd. Therefore, perturbative diagram resummation can not 
lead to Bjorken hydrodynamics: diagrams lead to a free streaming final state, while Bjorken 
hydrodynamics in the absence of confinement is an eternally interacting state. 

Of course, for realistic finite nuclei, when the distance between the nuclei becomes compara- 
ble to their radii, the expansion becomes three-dimensional: in this sense Bjorken hydrodynam- 
ics is not really an "eternally" interacting state. The onset of 3D expansion would indicate a 
transition from Bjorken to Landau hydrodynamics and is outside of the scope of this discussion. 
It would be very unlikely though for thermalization to happen at this late stage of the collision 
when the system becomes very dilute and the energy density starts falling much faster with 
time than it was in the Bjorken hydrodynamics stage of the collision. 

One may suggest that Feynman diagrams may lead to Bjorken hydrodynamics in some 
intermediate state preceding free streaming. Again this is not possible, since, as was shown 
above, Bjorken hydrodynamics does not lead to free streaming, and if the system enters a 
Bjorken hydrodynamics state it would never leave it. Therefore Bjorken hydrodynamics can 
not be an intermediate state of the system. 

Introducing massive quarks in the argument is not going to change the conclusions: if the 
quark mass is much smaller than the temperature, m q <C T, then it could be neglected yielding 
o ~ a 2 s /ct s T 2 ~ 1/T 2 and p ~ T 3 leading to the mean free path of Eq. (|13p. If quark mass is 
much bigger than the temperature, m q 3> T, then the particle density would still be dominated 
by gluons yielding p ~ T 3 and the interaction cross sections would also be dominated by gluon 
exchanges, giving cr ~ 1/T 2 , and we would recover Eq. (fT3~j) again. 

The argument presented above shows that Bjorken hydrodynamics is fundamentally dif- 
ferent, for instance, from the case of Landau hydrodynamics. In the latter case decoupling 
is possible even in the perturbative case. For simplicity let us consider an expanding three- 
dimensional spherical droplet with the energy density within the droplet independent of spacial 
coordinates and with the ideal gas equation of state (jBJ). Then for such a system 



Following the above steps we obtain, using Eq. pip, that T ~ 1/t and the mean free path is 

Xpert ~ — ~ (22) 

Since the radius of an ultrarelativistic droplet is R ~ t we see that \ per t ~ R and the sys- 
tem may decouple at late times. Therefore, thermalization in this system, leading to Landau 
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hydrodynamics description of its evolution, appears to be possible to describe using Feynman 
diagrams as an intermediate state between the creation of the system and eventual decoupling. 
Indeed to make a definitive conclusion on the subject of perturbative thermalization in 3D ex- 
panding system one needs to perform a careful diagrammatic analysis similar to the one carried 
out in pQ for the case of Bjorken hydrodynamics. 

One may worry that, for instance, Boltzmann equation has a diagrammatic interpretation 
[86J, i.e., it does sum up certain classes of diagrams in its collision term, and, nevertheless, it 
does lead to perturbative thermalization. However we would like to point out that Boltzmann 
equation leads to thermalization for a system of finite volume, such as a gas of particles in a 
box. If a system is not confined to a fixed volume, the question of thermalization of the system 
in the weak coupling limit due to Boltzmann equation dynamics is still an open question which 
has to be addressed on individual basis for each system. 

For the case of a system in a box, thermalization time in the weak coupling limit is compa- 
rable to the size of the box. In the diagrammatic language this implies that if we iterate the 
diagrams from the collision term of the Boltzmann equation to obtain a full resummed diagram 
describing the system, such diagram, covering the evolution of the system over a very long time 
scale in a fixed finite volume, should also include information on interactions of particles with 
the box in which they are confined, and, because of that, would not have an asymptotic "free 
streaming" final state. Therefore, it would not be a diagram in the vacuum and the above 
argument would not apply to it: one can think of the box as of an external field, which is 
present at all times preventing decoupling of particle interactions and not allowing them to free 
stream to spacial infinity. 

To conclude this Section let us reiterate our argument once more. Bjorken hydrodynamics 
j2j, in the weak coupling limit (in the absence of confinement effects), does not decouple: gluons 
and quarks continue interacting indefinitely, as long as one-dimensional expansion continues. 
Such system does not have a free final state, and, therefore can not be obtained from Feynman 
diagrams in vacuum, which, at any finite order in the coupling, and, therefore, to all orders as 
well, always do have a non-interacting final state. (Indeed Feynman diagrams in an external 
field, such as the graphs describing a system of particles in a box, do not have an asymptotically 
free final state: however, since, of course, there is no external field in heavy ion collisions, such 
diagrams would be irrelevant here.) 

2.2 Argument II 

Now let us assume that the system of quarks and gluons produced in a heavy ion collision reaches 
equilibration at some proper time r t h- Let us pick some proper time To after thermalization, 
r o > Ttht and set the QCD coupling constant g to be equal to zero for all times after r , i.e., 
g = for r > r . Below we will compare the g — > limits of the full theory (QCD) and of 
Bjorken hydrodynamics. 

First of all, from QCD standpoint, if the coupling goes to zero, g — > 0, the system would 
stop interacting and quarks and gluons would "free-stream" away. The energy density of a 
free-streaming quarks and gluons scales as 

e ~ - (23) 
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simply due to energy conservation in a one-dimensionally expanding system. (To obtain a 
different energy scaling with proper time, such as in Bjorken hydrodynamics case in Eq. (JHJ), 
the system needs to perform work, which obviously does not happen with a free-streaming 
system.) 

To analyze the g —>■ limit of Bjorken hydrodynamics we first note that it requires only 
conservation of the energy-momentum tensor (Eqs. (0) and Q) with an equation of state 
relating e and p. Since, in QCD, perturbative equation of state is known and is given by |H7] 



^ = 1 + ^ 2 ), (24) 

we naturally conclude that in g — > limit it reduces to the ideal gas equation of state (JHJ). 
Using that equation of state in Eq. (J3j) governing Bjorken hydrodynamics [2] leads to 

(25) 



r 4/3' 

Thus we have established that the g — > limits of a field theory, such as QCD, and of Bjorken 
hydrodynamics are different, given by Eqs. and (|2*5|l correspondingly. Therefore, pertur- 
bative QCD can never lead to ideal Bjorken hydrodynamics. 

The origin of the difference between the two limits is clear: the ideal gas equation of state 
(jSJ describes a gas of particles which do not interact with each other, but which interact with 
some thermal bath. This thermal bath could be a box the particles are confined to or some 
external field. Indeed there is no external thermal bath in the case of a heavy ion collision: 
hence Bjorken hydrodynamic approximation has a different g — > limit from the full theory. 

One may worry that putting the coupling constant to zero at proper time r > t would 
violate gauge invariance. Indeed this is true in general making our argument not quite rigorous. 
However, as was argued convincingly in the existing thermalization scenarios EEl ETJ EE 

EH1 EH]; by the time of possible thermalization (and even much earlier) the density of 
particles would become sufficiently low for the quarks and gluons to go on mass shell between 
the interactions. This would insure gauge invariance of the quark and gluon distributions at 
any given late proper time. 

To reconcile the g — > limit of Bjorken hydrodynamics from Eq. (|2*5j) with the more 
physically intuitive scaling of Eq. (|23|1 one has to include non-equilibrium viscous corrections 
[88, ESI- Introducing shear viscosity dependence into the ideal Bjorken hydrodynamics energy- 
momentum tensor would make it depend on three variables — e, p and rj, and is accomplished 
by adding a viscosity-dependent term to Eq. (j2J) . The modified energy- momentum tensor for a 
baryon-free fluid is [HH1 E2 IBS] 

T fW = {e + p)u ll u u -p g^ + rj {v ^u v + V„u M - jj A M „ V p u p ^j , (26) 

where V M = (g^ —u^Uy) d u and = g^ v —u^u u . The shear viscosity in perturbation theory 
for QCD with two light flavors is given by [72] 

rj = 86.473 ^ (27) 
9 ln(l/#) 
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and is, therefore, divergent in the g — > limit. Indeed this indicates importance of non- 
equilibrium corrections to Bjorken hydrodynamics, which was pointed out as a possibility al- 
ready in the original paper 2]. The non-equilibrium corrections become infinitely large in the 
g — ► limit implying a possible breakdown of Bjorken hydrodynamics. 

When viscous corrections are added to Bjorken hydrodynamics, the energy momentum 
tensor from Eq. (J2o^) can be written as (in the local fluid rest frame) 



( 



e(r) 






p{ T , 





T 3 r 









P(r) + 




2 2 

3 T 










(28) 



As follows from Eq. (J28)l . including viscosity-dependent additive term in the definition of T^ v 
reduces the longitudinal pressure. This may lead to the energy density scaling as shown in 
Eq. (J2EJ), since the pressure p on the right hand side of Eq. fl5J) is, in fact, the longitudinal 
pressure component. In other words, substituting Eq. into Eq. ((TJ) modifies Eq. (JSJ) to give 



e+[p-(4/3)»7/r] 



(29) 



One can see that shear viscosity tends to reduce the expression in square brackets on the right 
hand side of Eq. (|29~j) . making the scaling of energy density e with r closer to that of Eq. (J2BJ). 
Indeed if shear viscosity is very large this expression in square brackets may become negative: 
however, large viscosity needed for that implies that other non-equilibrium corrections to the 
energy-momentum tensor of Eq. (higher order derivatives of fluid velocity than the ones 
included in Eq. (J2EJ)) would also become important invalidating Eqs. (J2Hj) and (|29~|l . When 
departure from equilibrium becomes this large, it is not clear whether hydrodynamic approach 
would remain the best approximation to the relevant physics. Nevertheless, we have shown that 
non-equilibrium corrections tend to bring Bjorken hydrodynamics closer to the correct g — > 
limit of the full theory, indicating that the system at weak coupling is not in equilibrium. 

Here we would like to add a word of caution. As was shown in [T|, the most general energy- 
momentum tensor for a central collision of two ultrarelativistic nuclei of infinite transverse 
extent is given by 





rpfjLV 






V o 



P2[T) 









Pa(r) 




\ 







Pb(t)J 



(30) 



at z — 0. Therefore the most general T^ v depends only on three variables: energy density e, 
and transverse (p?) and longitudinal (^3) pressure components. The form of the tensors in Eqs. 
(|28|l and (J30|) is very similar: in fact, one can always express P2 and p% in terms of p and rj and 
vice versa. Equal number of degrees of freedom in viscous hydrodynamics of Eq. and in 
the most general energy-momentum tensor for the given collision geometry in Eq. (|3Tlj) makes 
it harder to distinguish the system really obeying the rules of viscous hydrodynamics from the 
system in some non-equilibrium state. 

Does our conclusion above indicate that ideal Bjorken hydrodynamics is impossible to 
achieve in heavy ion collisions? Possibly yes. However, one could also imagine the situa- 
tion where the thermalization and hydrodynamic behavior of the quark-gluon system are due 
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to non-perturbative effects, which disappear in the g — ► limit. In that case hydro dynamic 
description would only be correct in the strong coupling limit of the theory. An example of 
possible scaling of energy density could be given by the following ansatz 



r l+A cxp(-c/a s ) ' V > 

with A > and c > some non-perturbative constants. Indeed in the a s — ► limit Eq. (p!Tj) 
leads to Eq. (J2BJ)- It would also explain the absence of perturbative a s lnr corrections to the 
scaling of Eq. 1)23)1 observed in pQ : the function in Eq. 1)31)1 is non-analytic in a s and its Taylor 
series has zero coefficients for all non-zero powers of a s . The strong coupling limit of Eq. pip. 
a s — > oo, would give Eq. © thus recovering Bjorken hydrodynamics |2]. 

To conclude this Section let us summarize our argument presented here one more time: the 
difference of g — > limits of perturbative QCD and Bjorken hydrodynamics, given by Eqs. 
(|2H)l and correspondingly, indicates that it is impossible for perturbative evolution of the 
system to lead to Bjorken hydrodynamics. Non-perturbative ansatze of the type shown in 
Eq. (|3T)l are still possible though, indicating that success of hydrodynamic description of heavy 
ion collisions may truly be due to strong coupling effects, as appears to be demonstrated by 
RHIC phenomenology jHIHUHH] and as suggested in some thermalization models f90j. 



3 Non- Perturbative Contributions to Jet Quenching and 
Thermalization 

Below we will try to investigate the role of non-perturbative effects in heavy ion collisions by 
modeling them with an instanton ensemble. Instanton models have been very successful in 
describing many features of QCD vacuum. At the same time we are fully aware that instanton 
models do not capture all the properties of the QCD vacuum: the best-known example is the 
inability of instantons to explain confinement. Therefore, our discussion below describes only 
the contribution of instanton fields to the dynamics of heavy ion collisions, and does not include 
other non-perturbative (possibly strong coupling) effects, which may also be very important. 



3.1 Instanton Density in Color Glass Background 

The early stages of a heavy ion collision with proper times t<1/Q s , are dominated by strong 
(mostly classical) gluonic fields produced by the Color Glass Condensates in the colliding nuclei 
EH E2HH3J 12H 122 123 121 I22| ■ As was shown in [HOI, these strong gluon fields modify the vacuum 
instanton density [77) [TS] 



n {p) = 

to give (see also |79j ) 

n sat {p) = n {p) exp 



0.466e" L679jVe 1 / 2tt ' 
{N c -l)\{N c -2)\? \Mp), 



2N C 



(32) 



^p 4 



a.(JV? - 1) 



(A\ Gl v {x)Gl v {x) ~ G%{x)Gl v {x) \A) 



(33) 
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Here p is the instanton size, x is the instanton's position in space-time, G" is the field strength 
tensor of the produced gluon field and (A\ . . . \A) denotes averaging over nuclear wave functions. 
Eq. (j33| shows that instanton density may be significantly affected by the strong gluonic fields 
produced in the collision. 

The dominant gluon field produced in a heavy ion collision in the saturation/CGC frame- 
work is classical |2*T) 1221 123 I2H 123- Since the complete classical gluon field is known only 
numerically [21], following jHO] we will use the lowest order gluon field to calculate the corre- 
lators in Eq. (|3*5]). As was shown in [SU], the lowest order (o(g 3 )) gluon field produced in a 
nuclear collision, which was calculated in [22J, gives, in the forward light cone 

AiA 2 C F r d 2 k 



2^2 



(2tt) 2 S x S 2 tt 3 J (k 

dk+dk-dk'dk' e -^+*—^-*+-ik' + *— ik '- x + 

(2vr) 2 (fc 2 + iek ){k' 2 + iek' ) ' 1 } 

where A x and A 2 are atomic numbers of two nuclei which, for simplicity, are taken to be 
cylinders with their axes aligned with the collision axis and with cross-sectional areas S\ and 
5*2 correspondingly. In arriving at Eq. (|34|) one has to use the following formula (see Appendix 
B of JT], n, m are integers) 

dk + dk jll" k u X+ k\k™ = -2 7T 2 f4^r Jrn-n(k T T) (35) 

k + iek Q \2x- ) \ 2x + ) 

to replace each product k + k_ and k' + k'_ by k\j1 in anticipation of the integration over 
k + , k_, k' + , k'_. (If n = m the power of (k + kJ) n on the left hand side of Eq. (|3*3j) only brings in 
a power of (k^/2) n on its right hand side.) Here kx = \k\ and in Eq. k = —k [HHj. 

With the help of Eq. (|35J) we can perform the integration over k + , k_, k' + , k'_ in Eq. 
obtaining 

G%(x)G; u (x)) lo = -16 a * a C F ^± J ^4 {[Hk T r)f - [Mkrr)} 2 } . (36) 

Using the fact that, in the same lowest order approximation, the multiplicity of the produced 
gluons with transverse momentum kx, rapidity y and located at impact parameter b is |21 [ 12*2 } 

dN = 8a 3 s C F A, A 2 1_ k? 
d 2 kdyd 2 b tt S l S 2 k\ n A ' [ ' 

and neglecting the logarithm, we rewrite Eq. 1)36)1 approximately as 



[G%(x)G« u ( X )) lo * -2vr / **^JL_# {[Uk T r)f - [Uk T r)f} . (38) 

At the same time, as was shown in jSOj 

[G%(x)GUx)) t „ = 0. (39) 
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Combining Eqs. (|58j) and (j3Uj) with Eq. we write 

n sat {p) = n (p) eX p|--^- Ty 2 7 r/ {^Mf-[JiMf}}- 

(40) 

To analyze Eq. (|4()jl we will use the fact that the gluon transverse momentum spectrum 
due to classical fields scales as ~ l/h? for kT>Q s , as shown in Eq. (jTfjl . and is approximately 
constant (up to logarithms) for k T <Q s [2H 1221 123 I2H 1213 ■ This implies that the fc T -integral 
in Eq. (j4*Uj) is (approximately) dominated by k T w Q s . At early proper times, for r<l/Q s , the 
Bessel function J\{htr) ~ kxr/2 can be neglected compared to Jo(£>tt) ps 1 in Eq. (}4*0|) leading 
to instanton suppression as was pointed out in 

On the other hand, at late times, when t>1/Q s , we use the large-argument asymptotics of 
Bessel function to write 

[J (k T r)] 2 - [Ji{k T r)f « — — cos (Wr - fe T r > 1. (41) 

At these late times the exponent of Eq. (|4"0|) starts oscillating around with a decreasing 
amplitude. This means that instanton suppression by produced gluon fields disappears for 
r>l/Q s . 

This disappearance of suppression has a straightforward physical meaning. As was shown 
in pQ, any (classical and/or quantum) gluon fields in ultrarelativistic heavy ion collisions lead 
to the dominant contribution to the resulting energy density scaling as 

At proper times t>1/Q s , Eq. (j4*2^1 leads to the energy scaling of Eq. (J2EJ), corresponding to free 
streaming pQ. This means that by the proper time r> 1/Q S the gluon fields would linearize and 
propagate as plane waves. Plane wave gluons have equal chromo-electric and chromo-magnetic 
fields, E a = B a . Remembering that G a l = 2 (B a2 — E a2 ) we conclude that G a l w Oat times 
r> 1/Qs and suppression of instanton density disappears from Eq. (JHBJ)- 



3.2 Instanton Contribution to Jet Quenching 

Let us summarize the emerging picture of the collision. At early times, t<1/Q s , strong gluon 
fields are produced and instantons are suppressed. At later, but still comparatively early times, 
t>1/Q s , the gluon fields linearize and vacuum effects due to instantons are not suppressed 
anymore. The energy density of the linearized gluon fields is still quite high at r ~ 1/Q S , since, 
for r > 1/Q S 

o 3 

e ~ (43) 

as follows from Eq. (|42|) . (The factor of l/a s in Eq. (}4*3*|) is due to the dominance of classical 
gluon fields at the early stages of the collision [2U 1221 123 12111231221 123 •) From here on we will 
only perform parametric estimates to demonstrate the principle points: detailed calculations 
will be left for future analyses. 
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We assume that hadronization does not set in until the energy density of the gluon medium 
becomes of the order of ~ ^qcd- This energy density is obtained by requiring that we 
have one hadron with mass ~ ^qcd per unit volume in the rest frame of the hadron, which is 
roughly ~ 1/Aq CD . Equating the energy density from Eq. (f4*3"j) to Aq CD yields 

r h ^— . (44) 

Eq. (jH| assumes that one- dimensional expansion lasts indefinitely long. Of course in a realistic 
heavy ion collision one dimensional expansion lasts only up to r ~ R with R the nuclear 
radius. After that the expansion becomes three-dimensional and the hadronization energy 
density ~ Aq CD is achieved faster: we will return to this below. 

As we saw in Sect. 3.1, at early times, t<1/Q s , the non-perturbative instanton effects 
are suppressed and the vacuum is perturbative. At later times, t>1/Q s , instantons are not 
suppressed anymore and the vacuum is back to being non-perturbative. Deconfined gluons 
and quarks propagating through this vacuum would interact with the gluon (and quark) fields 
of the instantons. This may lead to energy loss by these partons in the instanton ensemble. 
Quarks and gluons, which, in a purely perturbative scenario, would have been propagating 
as plane waves for proper times 1/<3s<t < r^, start interacting with the non-perturbative 
vacuum fields loosing their energy. Again this phenomenon is much more pronounced in AA 
collisions, since in this case gluons and quarks spend a much longer time in the deconfined 
state interacting with the vacuum fields than they would in pp, pA or e + e~ collisions. Indeed 
after hadronization, when the gluons and quarks would become confined withing hadrons, the 
energy loss in vacuum fields would stop: hadrons are eigenstates of the QCD hamiltonian and 
they do not loose energy in propagation. 

The idea of partonic energy loss in the strong non-perturbative fields has been previously 
put forward by Shuryak and Zahed in jHIj. However, in their analysis the non-perturbative 
fields (sphalerons) were generated in a nuclear collision (see also [§2]), an d were not the QCD 
vacuum fields. 

To estimate the energy loss in the instanton vacuum fields by a propagating parton we will 
employ the results of calculations carried out in [HI EEHJ El CHI ECU EDI IBS ES] • Following [H2! 
we write the following expression for the energy loss of a quark created in a medium and having 
traveled distance z in that medium 

dE a q N r 1 , . r , Q Q da 



zp{z) J d 2 qq 2 T -^- : (45) 



dz 2 2 - a rK ' J ^ J d 2 q 

where p(z) is the density of the medium at point z. (Gluon energy loss is different from Eq. (J43J) 
only by a Casimir.) Here a is a constant depending on the rate of the medium expansion with 
time, < a < 1 j82j . Since we are interested in parametric estimates only, the exact value of a 
is not important to us. For the same reason the exact form of the (integrated) scattering cross 
section on a color charge in the medium in Eq. ()45|) is not important: we will just use the fact 
that it is 

"rfWj^-c* (46) 

Using such parametric estimates, we conclude that the integrated energy loss of a quark or 
gluon produced in a medium at proper time tq which travels through the medium up to some 
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hadronization time is 



- AE ~ a s / drrpir). (47) 



TO 



Below we will use Eq. (14 7|) to estimate the energy loss of a hard parton in the instanton vacuum 
fields and in the produced medium. 

Of course the energy density of the instanton vacuum is much lower (at early times) than 
the energy density of the produced medium. Therefore, one might assume the energy loss in 
vacuum fields to be negligibly small compared to the energy loss in the produced medium. 
However, for a wide range of functions p(r), the r-integral in Eq. (}4Tj) is dominated by late 
proper times r due to an extra factor of r in the integrand. Therefore, for the integrated energy 
loss, the density of the medium at late times is more important than the density of the medium 
at earlier times. Since, at late times, the density of produced medium becomes comparable 
to the instanton density in vacuum, the contributions of both the medium and the vacuum to 
partonic energy loss are quite comparable, as we will see below. 

The instanton density in the four- dimensional space-time is approximately [THj 

n « lfm" 4 « A 4 QCD . (48) 

Assuming that each instanton (or anti-instanton) provides at least one scattering center for the 
energy loss calculations of the type carried out in [1^1 HEl El CHI CHI H3 IH2 IHB] , Eq. gBJ allows 
us to put a lower bound on the resulting density of scattering centers at 

Pmst (r) « A 3 QCD . (49) 

Substituting Eq. pHj) into Eq. (f4Tj) yields the following energy loss in the instanton vacuum 
fields (r h > r ) 

-AE mst ~ alk% CD r 2 h . (50) 

To see whether the energy loss given by our estimates in Eq. (|5U() is large or small we want 
to compare it to the energy loss in the quark-gluon plasma. Indeed, as we argued above, this 
thermalized medium (QGP) can not be generated perturbatively. Below we will investigate 
whether parton-instanton interactions can lead to creation of a thermalized medium at early 
times. In the meantime we would just assume that a thermal medium with the energy density 
scaling given by Eq. © has been created in the collision. For this medium, requiring that the 
energy density should map onto Eq. (J4*3| at r = 1/Q S , we write 



Qi 

as (Q s r)V3- 



~ I TT^ZE - (51) 



Pq9p( T ) 3/4 1 ' ( 52 ) 



Remembering that if e ~ T 4 then p ~ T 3 we deduce from Eq. (151(1 the density of the produced 
medium 

«s /4 Q s T 

Substituting Eq. (|H2() into Eq. (|37|) gives us a parametric estimate of energy loss for a parton 
in quark-gluon plasma 

- AE qgp ~ a 9 s /4 Q 2 s T h . (53) 
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Now we can compare the energy loss in the instanton vacuum from Eq. (|5U|) to the energy 
loss in QGP from Eq. ([53)1 . First, assuming that nuclei have infinite transverse extent, such 
that the hadronization time is given by Eq. (JUJ, we substitute it into Eqs. (}50j) and (J53|) to 
obtain 

- AE inst ~ « s -^2- (54) 

A QCD 

and 

- A£ 99P ~ c^ 4 -^L (55) 

Since Q s 3> ^qcd and a s < 1 we conclude that 

|A£ inst | > |A£ g J, (56) 

i.e., that the energy loss in the instanton vacuum is actually somewhat greater than the energy 
loss in the produced medium! Since the saturation scale grows with energy and atomic number, 
Q 2 S ~ A 1 / 3 s A / 2 , the instanton-induced energy loss would increase with energy and the system's 
size faster than the energy loss in a medium. 

However, the conclusion of Eq. (JSljj) was reached for the one-dimensional expansion only. 
For realistic nuclei one-dimensional expansion stops at r ~ R. Since the hadronization time 
Th from Eq. (}4*4"|) is parametrically much larger than R, the one-dimensional expansion stops 
before the time is achieved. After that the expansion becomes three-dimensional. Assuming 
that free streaming continues in the 3D expansion, the energy density for the proper times after 
r m R is given by 

e = e(r = R) (^f (57) 
where e(r = R) could be obtained by putting r = R in Eq. (|4*Hjl . Eq. (jBTj) becomes 

Ql (58) 



a.R \t 



Equating the energy density from Eq. ()58|) to ~ Aq CD we obtain the hadronization time for 
a realistic nuclear collision 

rl D ~ R-^ - (59) 

If we use the hadronization time from Eq. (|59|) in Eqs. (|5()|) and (J53j) we get 

- AE inst ~ a 7 /' -5L (i2 Aqcd )4/3 (60) 

and 

-A^ p ~ a f/ 12 -^(/?A QOT ) 2 / 3 . (61) 

To compare Eqs. (f6T)|) and (jHTj) we first analyze them in the quasi-classical approximation: this 
regime seems to be adequate for mid-rapidity AA collisions at RHIC. As in the quasi-classical 
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McLerran-Venugopalan model Q 2 ~ A l ^h^ CD and since RAqcd ~ A 1 / 3 , we can conclude, 
neglecting the factors of ot s in the front, that for Eqs. (JoTIj) and (pTTj) 

\AE inst \ ~ A 7/9 > \AE qgp \ ~ A 13 / 18 . (62) 

Again, the instanton-induced energy loss appears to be slightly higher than the one in the 
plasma. Indeed, when quantum evolution corrections are included the saturation scale gets 
extra enhancement as Q 2 ~ s x ^ 2 : in that case \AE ggp \ from Eq. ()61|) would grow faster with 
energy than \AEi ns t\ from Eq. (|6()jl . eventually becoming larger. Still we have shown that the 
energy loss in the instanton vacuum is parametrically comparable to the energy loss in QGP even 
for realistic nuclear collisions, except for very high energies when Q s ^> Aqcd (RAqcd) 2 ^ 3 - The 
instanton energy loss from both Eqs. ()54|) and (J6*U|) also increases with energy since Q 2 S ~ s A//2 , 
which is in qualitative agreement with the absence of suppression of hadron production at lower 
(SPS) energy AA collisions and with the onset of suppression at higher (RHIC) energies. 

Indeed there is no real physical competition between the energy loss in QGP and in the in- 
stanton vacuum. If QGP is formed, instantons would be suppressed [93J, much like suppression 
in CGC that was discussed in Sect. 3.1. Here we calculate energy loss in QGP just to have a 
reference to compare instanton energy loss to. We need to do this since all our estimates in this 
Section are parametric and can not be used to produce numbers to compare to experimental 
data. 

Let us also note that a significant energy loss due to instanton fields only takes place in AA 
collisions. For instance, in pA collisions energy loss is given by Eqs. (j6*Uj) and ()61|) with the 
nuclear radius R replaced by the proton radius, r£ ~ r p ~ 1/Aq C d and the nuclear saturation 
scale Q s replaced by the proton saturation scale, which, in the quasi-classical approximation is 
roughly ~ Aqcd- After performing these substitutions we get 

- AEi? st ~ Aqcd ~ ~AE^ p , (63) 
which is very small, much smaller than the energy loss in AA collisions given by Eqs. (}60|) and 

flHU). 

To summarize, we have succeeded in showing that the energy loss in the vacuum fields 
modeled by an instanton ensemble is appreciably high, being parametrically comparable to 
the energy loss expected in hot quark-gluon plasma. Once again, let us reiterate that such 
conclusion is only possible due to a peculiar form of the integrated energy loss in Eq. (jUJ), 
which was derived in |82[ I83j. making AE mostly sensitive to the color charge density at late 
proper times. The obtained instanton energy loss increases with both center of mass energy of 
the collisions and with the size of the system. 

3.3 Instanton Contribution to Thermalization 

Finally, let us address the following question: if the interactions of the quarks and gluons with 
the instanton vacuum fields are sufficiently strong to generate a significant energy loss, as shown 
in Eqs. (J54J) and ()60|). should not the same interactions lead to thermalization of the quark- 
gluon medium? This is a difficult question which, in general, requires understanding to what 
extent interactions of partons with the instanton fields would trigger interactions of partons 
with each other, which, in a purely perturbative approach PHQ], would stop by r « l/Q s . 
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Here we will postpone this important question for further investigation, which is beyond the 
scope of this work, and concentrate on parton-instanton interactions only. 

As was shown in j2H [l] the perturbative gluon fields are produced by r « 1/Q S with zero 
longitudinal pressure component. Namely, at z = and at r « 1/Q S the energy- momentum 
tensor looks (approximately) like T^ u = diag(e,p,p, 0). Thermalization of the system must 
be accompanied by generation of the non-zero longitudinal pressure component comparable to 
the transverse pressure components jHS]- This is a necessary, though not sufficient, condition 
for thermalization. Let us estimate when such a longitudinal pressure component would be 
developed in parton-instanton interactions: that would give us a lower bound on thermalization 
time. 

Avoiding all the subtleties of parton-instanton interactions jSI] we assume that the upper 
bound on the cross section of a parton-instanton interaction is given by the instanton size 
squared, a inst ~ (p) 2 ~ 1/Aq CD [TS]. The mean free path of a parton in the instanton 
medium with the density given by Eq. (f4T)|) is then 

1 

Xinst = 1/Aqod. (64) 

Pinst &inst 

In each parton-instanton interaction, the momentum transfer to the parton would be of the 
order or 1/ (p) ~ Aqcd- After N rescatterings the parton's longitudinal momentum gets broad- 
ened to approximately 

Pz ~ VN Aqcd- (65) 

For the system to become isotropic we demand that p z ~ pr ~ Q s , which gives N ~ Q 2 S / A 2 q CD . 
The time it takes to have N ~ QII^qcd parton-instanton rescatterings is roughly given by 

O 2 

Tth ~ KnstN ~ -f*-. (66) 
A QCD 

This is a very late thermalization time. Even in McLerran-Venugopalan model, where Q 2 ~ 
A 1 / 3 Aq CD , this thermalization time becomes comparable to the nuclear radius, Tth ~ which 
corresponds to the end of the one-dimensional expansion. Certainly such thermalization time 
is too late for hydrodynamic models describing RHIC data [H E] • Including energy dependence 
in Q s would only make Tth larger, leading to thermalization time which increases with energy. 
Indeed a numerical prefactor in Eq. (|66p. when calculated, may change the actual numerical 
results somewhat. Nevertheless we have demonstrated that, at least for very large nuclei, 
instanton vacuum can significantly contribute to jet quenching, as was shown in Eq. ()6()jl . but 
is not likely to bring in fast thermalization. 



4 Conclusions 

Above we have presented two arguments demonstrating that perturbative thermalization lead- 
ing to ideal Bjorken hydrodynamics is impossible in heavy ion collisions. The only possible way 
to reconcile this with the phenomenological evidence for fast equilibration of the quark-gluon 
system produced in the collision |U E] is to suggest that thermalization takes place through 
some non-perturbative processes. It might be reasonable to assume that non-perturbative ef- 
fects become important at proper times ta ~ 1/Aqcd in ^4^4 collisions, as they are known 
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to do in pp, pA and e + e~ collisions. Non-perturbative effects involve QCD strings jHl] with 
the string tension a pa 1 GeV/fm. Gluons and quarks produced in the initial stages of the 
collision carry a typical transverse momentum of the order of Q s , which is Q s = 1 4- 1.4 GeV at 
RHIC [51]. This number is comparable to <jt\ pa 1 GeV, indicating that QCD string effects 
may become important at proper time r\, introducing strong-coupling interactions between 
the partons and, possibly, bringing the system towards fast thermalization. The corresponding 
time scale ta pa 1/Aq CD pa 1 fm/c appears to be consistent with the thermalization time of 
approximately 0.5 fm/c required by hydrodynamical models [H |3] to describe RHIC elliptic 
flow data [Zj. 

Non-perturbative thermalization, either through instantons or through QCD string interac- 
tions, would not violate our argument of impossibility of perturbative thermalization presented 
above and in pQ. As we have argued in Sect. 2.2, thermalization in heavy ion collisions could 
only be a strong coupling phenomenon: such effects lie outside of the analysis carried out in 
p. The small coupling a s used in the energy loss discussion of Sect. 3.2 is due to parton 
momentum broadening in the non-perturbative instanton vacuum and does not infringe on the 
non-perturbative nature of the vacuum. Indeed, the thermalization scenario discussed in Sect. 
3.3 is essentially partonic, but it requires non-perturbative background instanton fields, which 
are also outside of the perturbative analysis of [T]. 

What appears to be unclear in non-perturbative thermalization scenarios with either instan- 
ton or string interactions is the energy dependence of thermalization time. It is not clear how 
string-induced interactions "know" about the center-of-mass energy of the collision. Instanton- 
mediated thermalization appears to have thermalization time (I66|) increasing with energy. This 
contradicts the fact that thermalization time required by phenomenology at RHIC is short, 
while it is not clear whether thermalization has ever occurred at the SPS implying a much 
longer thermalization time at lower energies. 

Non-perturbative effects appear to be needed to describe jet quenching in heavy ion colli- 
sions as well. As we discussed above in Sect. 3.1, in the purely perturbative scenario of the 
collision, by the proper time r ~ l/Q s the gluon fields would linearize and propagate as plane 
waves. (Similar results have been obtained for quark fields.) This result implies that, as long 
as perturbative interactions are concerned, nothing happens to the gluons after r ~ l/Q s - 
Therefore, the resulting gluon spectra would be completely determined by the (very) early 
stage dynamics. This would be in contradiction to the observation of Cronin enhancement at 
mid-rapidity in d + Au collisions at RHIC jHU 0U] and the observation of suppression of hadron 
production at mid-rapidity in Au + Au collisions jUS] : all the known perturbative production 
mechanisms imply that, since hadrons produced at mid-rapidity in d + Au collisions exhibit 
Cronin enhancement, so should the partons in the early stages of Au + Au collisions. Therefore, 
to reconcile RHIC data with the theory one has to assume that strong final state interactions 
take place in AA collisions, which could only be non-perturbative as was demonstrated by the 
above analysis. Such effects could be due to non-perturbative thermalization leading to forma- 
tion of QGP in which the partons would loose energy. Alternatively, as we have shown above, 
partons can loose energy in the non-perturbative vacuum fields, to which they are exposed until 
the hadronization time in heavy ion collisions. 

To summarize we note that we have demonstrated above that interactions in the final stages 
of heavy ion collisions can only be non-perturbative in nature. They may or may not lead to 
rapid thermalization, but they are likely to significantly contribute to jet quenching. 
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